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Precision  low  elevation  angle  radar  observations  of  orbiting  spheres  were  used 
to  evaluate  the  utility  of  surface  based  refraction  techniques.  Elevation  angle 
error  measurements  were  obtained  from  differences  between  the  precision  radar 
measurements  and*  the  true  or  absolute  position  of  the  orbiting  sphere  obtained 
from  a multiple  pass  orbit  determination  using  data  from  the  same  radar  system. 
The  results  of  the  refraction  effects  study  show  that  the  average  observed  ele- 
vation angle  error  caused  by  refraction  is  identical  with  the  average  calculated 
elevation  angle  error  obtained  from  ray  tracings  using  radiosonde  data  to  within 
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* the  measurement  accuracy  of  the  radar  system.  The  results  also  reveal  a 
relatively  large  uncorrectable  component  of  the  elevation  angle  error  caused  by 
small  scale  refractive  index  gradients.  This  uncorrectable  component  is  of  the 
order  of  30  mi 1 li degrees  (radeg)  at  a true  elevation  angle  of  1 degree  and  4 mdeg 
at  10  degrees. 
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EVALUATION 


1.  This  document  is  the  Final  Report  on  the  contract  which  over  the  period 
from  January  1977  to.  May invest! gated  tropospheric  effects  on  low 
elevation  angle  ray  paths.  The  report  describes  the  data  collection, 
verification  and  analysis  leading  to  a set  of  final  conclusions  which  are 
extremely  valuable  to  radar  system  designers  and  users.  Today's  precision 
radar  systems  require  accurate  tropospheric  refraction  corrections  to 
account  for  the  time  delay  and  ray  path  bending  due  to  the  lower  troposphere. 
The  question  arises:  Is  the  state  of  the  art  in  tropospheric  refraction 
predictions  such  that  further  advances  in  radar  equipment  accuracy  will  yield 
little  or  no  system  accuracy  improvement  due  to  tropospheric  prediction  errors? 
The  results  of  this  contract  provide  the  answer.  For  elevation  angles  below 
twenty  degrees,  the  fluctuations  due  to  large  and  small  scale  atmospheric 
processes  result  in  residual  errors  which  cannot  be  reduced  by  either  surface 
based  or  radiosonde  observations.  These  residual  errors  are  larger  than 
predicted  by  simulation  and  result  in  less  improvement  in  the  prediction  of 
tropospheric  refraction  errors  using  surface  data  than  expected. 

2.  The  inescapable  conclusion  of  this  effort  is  the  existence  of  a lower 
bound  in  the  prediction  accuracy  of  tropospheric  refractive  effects  which 
cannot  be  crossed  without  a vast  increase  in  the  knowledge  of  the  instant- 
aneous state  of  the  lower  atmosphere.  The  use  of  surface  refractivity  as  a 
predictor  results  in  errors  very  close  to  the  lower  bound  and  this  accounts 
for  the  relatively  small  improvement  derived  from  the  use  of  more  complex 
refraction  prediction  schemes. 

3.  The  results  of  this  contract  will  be  of  interest  to  radar  system  users 
and  designers  who  are  contemplating  systems  with  angular  accuracies  approach- 
ing the  lower  bound  described  above. 
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1 . INTRODUCTION 


1.1  Objective 

The  objective  of  this  contract  is  an  evaluation  of  the  adequacy 
of  surface  based  methods  for  the  estimation  of  refraction  effects 
at  low  elevation  angles.  A number  of  simulations  of  refraction  effects 
have  been  performed  using  ray  tracing  techniques  which  have  resulted 
in  procedures  for  the  correction  of  refraction  induced  errors  (Bean 
et  al,  1 960;  Crane,  1976a).  Adequate  experimental  verification  of 
the  correction  procedures,  however,  was  not  available.  As  one  of  the 
end  products  of  a large  scale  radar  propagation  experiment,  Lincoln 
Laboratory,  MIT  amassed  a number  of  precision  radar  tracks  of  orbiting 
spheres  sufficient  for  the  testing  of  refraction  correction  algorithms 
(Evans,  1973a;  Crane,  1976b).*  These  data  were  analyzed  under  this 
contract  to  experimentally  determine  refraction  induced  elevation 
angle  errors  at  low  elevation  angles  and  to  evaluate  the  utility  of 
a statistically  derived  correction  procedure  based  on  the  use  of 
surface  refractivity  observation, 

1.2  Summary  of  Results 

Precision  radar  tracks  of  orbiting  calibration  spheres  were  analyzed 
to  provide  statistical  estimates  of  refraction  induced  elevation  angle 
errors  and  the  residual  errors  remaining  after  correction  based  on 
surface  refractivity  observations.  The  elevation  errors  were  estimated 
from  differences  between  radar  observed  sphere  positions  and  the  ’’true" 
or  absolute  positions  of  the  sphere  obtained  by  calculation  from  the 
orbital  elements  for  each  sphere.  Precise  reference  orbits  were  estab- 
lished for  each  sphere  using  multiple  pass  radar  observations  and  a 
specially  developed  non-real-time-precision  orbit  determination  program, 
ORBFIT  (see  section  IV,  Evans,  1973b;  TRW,  1966). 


* Radar  tracking  data  and  precise  orbit  data  were  provided  courtesy 
of  Dr.  J.V.  Evans  of  Lincoln  Laboratory. 
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The  results  of  this  refraction  analysis  study  are  presented  in 
Tables  1 and  2.  The  true  elevation  angle  is  the  straight  line  (unre- 
fracted) direction  to  the  sphere.  The  observed  elevation  error  is  the 
refraction  induced  pointing  angle  error;  the  radar  observed  the  sphere 
at  an  apparent  elevation  angle  too  high  by  the  observed  elevation  angle 
error.  The  expected  elevation  angle  error  is  the  result  of  computer 
simulation  by  ray  tracing.  The  average  observed  and  expected  values  dif- 
fer by  at  most  3.2  percent  and  are  generally  within  2 percent  for  eleva- 
tion angles  below  10  degrees.  The  results  show  excellent  agreement 
between  observation  and  simulation;  the  differences  are  within  the  mea- 
surement error  of  the  Millstone  L-band  radar  system.  Based  on  these 
results,  it  is  evident  that  surface  refraction  correction  procedures  can 
reduce  the  average  elevation  angle  error. 

Table  2 presents  the  root  mean  square  (rms)  deviations  of  the 
elevation  errors  about  t’e  mean  value  or  about  the  estimated  elevation 
angle  correction  calculated  using  the  simultaneously  observed  value 
of  surface  refractivity.  Earlier  work  by  Crane  (1976b)  showed  that 
short  term,  random  fluctuations  in  elevation  angle  were  of  the  same 
order  of  magnitude  as  the  expected  day-to-day  variations  (expected 
residual  after  correction) . In  that  work,  only  deviations  about 
the  detrended  elevation  angle  were  considered;  bias  or  average  refraction 
effects  were  ignored  and  only  the  random  fluctuations  were  analyzed. 

In  this  study,  the  rapidly  fluctuating  components  of  the  elevation  angle 
error  were  removed  by  averaging  and  only  the  slowly  varying  or  bias 
values  are  considered.  The  observed  deviations  are  from  one  20-second 
observation  interval  to  the  next  and  from  one  satellite  pass  to  the  next; 
they  represent  the  random  variation  of  the  bias  component  of  the  refrac- 
tion error.  The  observed  fluctuations  about  the  mean  range  from  equal 
to  nearly  1.4  times  the  expected  (simulated)  value,  as  shown  in  Table  3. 
The  observed  residual  after  correction  ranged  from  2.5  to  7.4  times  the 
expected  value;  the  reduction  in  rms  deviation  was  by  less  than  a factor 
of  1.7  while  the  expected  reduction  was  by  a factor  of  5.9  or  better  at 
elevation  angles  above  1°. 

The  surface  correction  procedure  reduced  the  rms  deviations  by  as 
much  as  was  possible  for  the  given  set  of  surface  refractivity  values. 

The  correlation  between  the  residual  elevation  angle  after  correction 
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and  the  surface  refractivity  value  used  for  correction  was  not  signifi- 
cantly different  from  zero  (see  Table  1;  the  confidence  intervals  are 
for  67  percent  included  or  1 standard  deviation,  1 a).  This  result, 
obtained  from  a sample  of  over  1S00  bias  elevation  angle  error  observa- 
tions, shows  that  although  the  correction  technique  did  well  on  average, 
considerable  error  may  still  exist  for  a particular  pointing  angle  obser- 
vation. The  surface  correction  scheme  did  not  work  as  well  as  expected 
based  on  the  simulation  analysis.  Since  the  earlier  analysis  of  Haystack 
observations  (Crane,  1976b)  indicated  that  the  random  fluctuations  about 
a detrended  elevation  angle  value  were  of  the  order  of  the  expected 
residual  values  and  that  increased  fluctuations  should  occur  when  larger 
observing  intervals  in  either  time  or  elevation  angle  were  used,  the 
larger  random  fluctuations  for  the  larger  elevation  angle  spans  of  the 
Millstone  observations  were  expected.  These  larger  fluctuations  are 
random  in  nature  and  not  accessible  to  correction  using  surface  values. 
The  observed  rms  deviations  after  correction  represent  the  limitations 
of  surface  correction  procedures. 

Subsequent  sections  of  this  report  provide  a brief  background  for 
the  measurement  program,  a description  of  the  observations  and  ray  trac- 
ing simulations  and  an  analysis  of  the  results.  Section  2 places  the 
measurement  program  in  context.  Section  3 describes  the  Millstone  Hill 
L-band  measurement  program.  Section  4 describes  the  simulation  proce- 
dure and  the  ray  tracing  procedure  used  for  the  calculation  of  refraction 
effects.  Section  5 presents  the  comparison  between  observations  and 
expectations  and  Section  6 summarizes  the  conclusions. 
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2.  BACKGROUND 


2.1  Prior  Studies 

Radar  systems  operating  within  the  troposphere  have  long 
experienced  measurement  errors  caused  by  spatial  and  temporal  variations 
in  the  index  of  refraction  of  the  air.  Although  the  phenomenon  was  well 
known,  refraction  induced  errors  were  not  of  consequence  for  most  radar 
systems  operating  prior  to  1960  with  the  exception  of  height  finders. 

A number  of  simulation  studies  were  made  in  the  late  1950' s to  evaluate 
the  seriousness  of  refraction  induced  errors  and  to  provide  possible 
procedures  for  refraction  correction  (Bean  and  Thayer,  1954,  Bauer  et  al, 
1958;  Bean  et  al,  1960).  From  these  studies  it  became  evident  that  the 
calculated  values  of  ray  bending  and  of  elevation  angle  error  were  highly 
correlated  with  the  surface  value  of  the  index  of  refraction  or  surface 
refractivity . Bean  et  al  [1960]  provided  tables  for  the  statistical 
correction  of  elevation  angle  errors  using  surface  refractivity  measure- 
ments. The  correction  values  were  estimated  using  a linear  relationship 
between  elevation  angle  error  (or  bending)  and  surface  refractivity.  The 
coefficients  for  use  in  these  relationships  were  derived  from  linear 
regression  analysis  of  simulated  errors  and  surface  refractivity  values. 

The  ray  tracings  used  for  the  simulations  were  calculated  based  on  a 77 
profile  sample  representing  six  different  refractivity  profile  types  for 
13  climate  regions. 

The  use  of  the  statistical  prediction  technique  had  not  been 
experimentally  verified.  Bean  and  Thayer  [1963]  reported  on  the  use  of 
the  statistical  model  and  provided  samples  of  comparisons  between  simula- 
tions and  observations.  Their  most  convincing  comparison  with  experiment 
was  with  the  measurements  performed  by  Anway  [1963]  using  a 2-cm  wave- 
length radio  sextant.  The  results  reported  by  Anway  differed  from  the  mean 
values  of  elevation  angle  error  estimated  using  the  77  profile  sample 
by  less  than  the  standard  deviation  of  the  observations.  Anway  also 
reported  on  the  correlation  between  observed  elevation  angle  error  and 
surface  refractivity.  He  calculated  a correlation  coefficient  of  0.87 
for  observations  at  an  apparent  elevation  angle  of  8°,  sufficient  to 
reduce  the  rms  deviations  by  a factor  of  2.8  using  the  surface  correction 
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procedure.  The  measured  congelation  coefficient  was  not  as  high  as  the 
value  calculated  by  Bean  et  al,  a value  in  excess  of  0.99  at  an  apparent 
elevation  angle  of  8°.  However,  the  results  were  very  encouraging  and 
refraction  correction  based  on  surface  refractivity  values  have  been  used 
in  radar  system  operation  at  low  elevation  angles.  Recently,  in  an 
Investigation  of  the  effects  of  large  scale  horizontal  variations  in 
refractive  index  profiles,  Gallop  and  Telford  (1975)  simulated  the  effects 
of  horizontal  inhomogeneity  and  reported  that  correction  based  on  only 
surface  refractivity  measurements  would  not  be  as  successful  as  expected 
for  a horizontally  homogeneous  atmosphere.  They  suggested  that  the 
effects  of  horizontal  variation  would  be  small  at  a 5*  elevation  angle 
and  negligible  at  10°,  not  sufficient  to  explain  the  Anway  results. 


2.2  Millstone  Hill  Radar  Propagation  Study 


In  1969,  Lincoln  Lauoratory  began  a joint  study  with  the  Bell 
Telephone  Laboratory  to  investigate  the  limitations  imposed  by  naturally 
occuring  propagation  effects  on  defensive  radar  systems  operating  at  UHF 
(Evans,  1973a).  The  initinl  objective  of  the  joint  study  was  to  evaluate 
the  effects  of  the  auroral  region  of  the  ionosphere  on  radar  system 
performance.  Refraction  effects  associated  with  the  auroral  region 
were  expected  to  be  important  at  UHF  but  not  at  L-band.  Accordingly,  a 
two  frequency  measurement  program  was  devised  to  use  the  Millstone  Hill 
L-band  radar  for  observations  of  UHF  beacon  satellites  to  provide 
reference  position  information  and  a newly  developed  UHF  tracker  on  the 
same  antenna  for  observations  of  the  angle-of-arrival  of  the  signals 
from  the  beacons.  As  the  program  developed  and  more  stringent  accuracy 
requirements  were  placed  on  locating  the  absolute  position  of  the  beacon, 
the  use  of  the  L-band  radar  for  reference  position  determination  based 
on  single  radar  observations  was  abandoned. 

An  extensive  development  program  was  undertaken  to  provide  a mea- 
surement and  analysis  system  capable  of  accurately  determining  the  orbits 
of  the  beacon  satellites.  With  precise  orbital  information,  the  abso- 
lute location  of  the  satellite  could  be  determined  and  tracking  pertur- 
bations induced  by  the  ionosphere  could  be  determined  by  comparing  the 
observed  position  with  the  calculated  position.  To  accomplish  the  task 
of  obtaining  an  accurate  reference  orbit  for  a satellite,  the  satellite 
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was  observed  over  many  successive  passes  over  the  radar  station  and  the 
orbital  equations  were  solved  using  a non-real-time-precision  orbit 
determination  program  to  establish  the  best  fit  reference  orbit  consis- 
tent with  orbital  mechanics  and  a large  number  of  radar  observations. 

An  extensive  radar  calibration  program  was  undertaken  to  provide  models 
for  all  known  sources  of  errors  in  the  L-band  radar  observations.  The 
errors  had  to  be  modeled  and  corrected  prior  to  use  in  the  orbit  fitting 
program  to  provide  the  best  possible  data  for  the  determination  of  the 
reference  orbit  (Evans,  19731). 

As  a part  of  the  L-band  radar  calibration  and  error  modeling  pro- 
gram, a surface  refractivity  correction  scheme  was  developed  that  was 
tailored  to  the  Millstone  radar  climate  (northeastern  Massachusetts). 

The  statistical  procedure  developed  by  Bean  et  al  (1960)  was  based  on  a 
limited  number  of  profiles  representing  extreme  conditions  in  a number 
of  different  climate  regions.  A new  regression  analysis  was  undertaken 
(Evans,  1969;  Crane,  1976a)  using  radiosonde  data  from  two  months 
(February  and  August)  of  each  of  three  years  for  Albany,  New  York,  a 
location  with  a climate  similar  to  Millstone  Hill.  The  new  statistical 
parameters  were  expected  to  be  better  than  the  earlier  model  because  the 
parameters  were  from  a statistical  analysis  using  a sample  population 
which  represented  the  population  to  be  expected  at  observation  time. 

The  surface  refractivity  correction  procedure  was  used  without  experi- 
mental verification  throughout  the  joint  study  measurement  program, 
1970-1973.  The  regression  coefficients  for  use  at  low  elevation  angles 
were  modified  in  1972  to  reduce  differences  between  the  statistically 
predicted  bending  values  and  ray-tracing  calculations  for  four  soundings 
obtained  with  a nearby  radiosonde. 

The  joint  propagation  study  was  concluded  in  March  1973  and  a fol- 
low-on tropospheric  refraction  effects  study  was  begun  by  Lincoln  Lab- 
oratory for  the  U.S.  Army  Advanced  Ballistic  Missile  Defense  Agency 
(ABMDA) . The  tropospheric  effects  measurement  program  was  conducted  in 
two  phases,  the  first  a series  of  L-band  radar  observations  conducted 
in  sixteen  tracking  sessions  between  1 October  1973  and  30  September  1974 
for  the  purposes  of  measuring  the  bias  refraction  errors  and  evaluating 
the  surface  refractivity  correction  procedure  and  the  second  a series  of 
measurements  of  the  random  fluctuations  in  amplitude  and  angle-of-arri val 
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caused  by  refractive  index  inhomogeneities  in  the  lower  troposphere. 

The  random  fluctuation  or  scintillation  observations  were  made  during 
1975  and  have  been  reported  by  Crane  (1976b).  The  L-band  observations 
were  processed  through  the  determination  of  the  reference  orbits  for 
each  of  the  spheres  tracked  during  the  last  four  tracking  sessions  but 
the  data  were  not  analyzed  further  due  to  changing  priorities  both  in 
the  Army  and  at  Lincoln  Laboratory. 

This  report  presents  the  results  of  the  analysis  of  the  bias  refrac- 
tion data  obtained  during  the  last  two  tracking  sessions,  during  Septem- 
ber 1974.  The  completion  of  the  analysis  has  been  supported  by  the  U.S. 
Air  Force,  Rome  Air  Development  Center  to  provide  experimental  verifica- 
tion of  the  use  of  surface  refraction  correction  techniques  for  the  cor- 
rection of  elevation  angle  errors  at  low  elevation  angles. 
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3.  PRECISION  REFRACTION  EFFECTS  MEASUREMENT 


3,1  L-band  Radar  Systems 

The  Millstone  Hill  L-band  radar  system  has  been  described 
by  Ghiloni  [1973],  The  radar  system  is  a high  powered,  monopulse 
tvacker.  For  satellite  tracking  operations,  the  radar  system  was 
nominally  operated  at  a 3.3.  MW  peak  power  and  100  KW  average  power. 

A 2-millisecond  pulse  was  used  at  n pulse  repetition  frequency  of  15  pps. 
The  radar  utilizes  a Cassigrain  antenna  system  with  a 25.6m  (84-ft.) 
diameter  aperture.  At  the  operating  frequency  of  1295  MHz,  the  antenna 
has  a gain  of  about  47  dB  with  a 0.7  degree  half  power  beamwidth  (one- 
way). The  transmitted  polarization  was  right-hand  circular. 

The  receiver  system  employs  a 12-horn  monopulse  tracking  feed 
to  develop  the  sum  and  difference  signals  required  for  detection  and 
tracking.  The  left-hand  circular  polarization  sum  channel  signals  were 
processed  through  low  noise  RF  amplifiers,  several  stages  of  IF  processing 
and  were  both  coherently  and  incoherently  detected.  The  in-phase  and 
quadrature  signal  amplitude  values  and  the  normalizing  amplifier  AGC 
output  were  digitally  recorded  for  use  in  computer-aided  real-time  sat- 
ellite tracking  (using  the  CAST  program,  see  Evans  1973b)  and  for  post 
mission  analysis.  The  elevation  and  azimuth  difference  signals  in-phase 
with  the  sum  signal  were  normalized  using  the  AGC  output  to  derive 
pointing  angle  error  values  which  were  digitized  and  input  to  the  on-site 
SDS  9300  computer  for  real-time  tracking  and  storage  for  post  test  analysis 
The  radar  signals  were  processed  in  real-time  and  used  to  develop 
tracking  commands  to  steer  the  antenna,  the  range  gate,  and  the  Doppler 
frequency  offset  to  control  one  of  the  IF  frequencies.  The  antenna  posi- 
tion encoder  values  were  also  processed  by  the  computer  and  used  to 
generate  the  new  pointing  commands;  the  tracking  loops  were  closed  through  1 
the  computer. 

Auxilliary  data  for  pointing  error  correction  were  also  sam- 
pled and  processed  through  the  computer.  The  surface  refractivity  was 
measured  by  a microwave  refractometer  mounted  on  the  outside  wall  of  the 
radar  equipment  shelter  attached  to  the  antenna  and  it  was  positioned  at 
roughly  the  height  of  the  phase  center  of  the  antenna  when  operating  at 
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low  elevation  angles.  The  refractivity  values  were  averaged  for  1000 
seconds  prior  to  recording  to  remove  small  scale  turbulent  fluctuations 
in  the  refractivity  values.  The  tilt  of  the  antenna  tower  was  also 
measured  and  recorded  for  use  in  removing  the  effects  of  mechanical 
deformation  of  the  antenna  supporting  structure.  The  deformation  of  the 
tower  caused  by  solar  heating  often  exceeded  20  mdeg  and  required  compen- 
sation to  provide  angle-of-arrival  measurements  of  sufficient  accuracy 
to  determine  refraction  effects.  Tower  tilt  was  measured  by  two  ortho- 
gonally oriented  electronic  level  sensors,  the  principal  tilt  measure- 
ment was  in  the  elevation  plane  of  the  antenna. 

The  Computer  Aided  Satellite  Tracking  program  (CAST)  was  an  inte- 
gral part  of  the  L-band  radar  tracking  system  used  for  the  low  elevation 
angle  refraction  measurements.  This  program,  resident  in  the  SDS  9300 
computer,  closed  the  tracking  loops  used  to  steer  the  antenna,  range 
gate,  and  Doppler  offset . Single  L-band  radar  returns  were  not  suffi- 
ciently noise-free  for  precision  satellite  tracking.  Accordingly,  a 
tracking  procedure  was  devised  to  use  the  radar  observations  to  correct 
the  calculated  satellite  ephemeris  values;  the  radar  data  were  smoothed 
and  used  to  differentially  adjust  the  predicted  satellite  position.  In 
this  way,  a significant  amount  of  smoothing  could  be  used  to  reduce 
radar  measurement  uncertainty.  The  CAST  program  also  removed  known, 
modeled  measurement  errors  caused  by  refraction,  tower  tilt,  encoder 
offsets,  feed  droop  (change  in  the  electrical  axis  of  the  antenna  rela- 
tive to  the  mechanical  axis  caused  by  mechanical  deformation  of  the 
antenna  as  the  antenna  elevation  angle  is  changed),  beam  collimation 
and  non-perpendicularity  of  the  elevation  and  azimuth  axes.  The  smoothed 
differences  between  the  ephemeris  values  and  the  best  estimate,  corrected 
radar  position  observations  were  used  to  correct  the  ephemeris  values  to 
provide  the  tracking  commands. 

The  raw  data  were  corrected  to  provide  best  estimate  azimuth  and 
elevation  angle  values.  The  corrections  used  were: 

AZ'  = AZ  + ABIAS  + ASLOPE  x AZ 

enc  r r enc 

AZ"  * AZ*  + sin(TILT)  x sin  (TILTAZ  - AZ»)  x tan  (EL  ) 

EL'  = EL  - sin (TILT)  x cos (TILTAZ  - AZ') 
cnc 

EL"  = EL'  + [1  - cos (TILT)]  x tan(ELenc) 
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EL"'  - EL"  - EBIAS  + AEL  - ELS LOPE  x EL" 

EL""  » EL'"  - (A  + B x N ) 

AA2  = [ATRAV  + COLL  - SKEW  x sin (EL" ') ] /cos  (EL'") 

AZ"'  » AZ"  - ABIAS  + AAZ 

where  AZ^„„  and  EL^„„  are  the  azimuth  and  elevation  angle  encoder  values, 
enc  enc 

ABIAS^  and  ASLOPE  provide  the  corrections  for  known  deficiencies  in  the 
r r 

azimuth  encoder  system,  TILT  and  TILTAZ  are  the  principal  and  orthogonal 
tilt  meter  l-eadings,  EBIAS  and  ELSLOPE  correct  for  deficiencies  in  the 
elevation  encoder  system  and  feed  droop,  Nf,  is  surface  refractivity, 

AEL  and  ATRAV  are  the  elevation  and  traverse  (azimuth)  error  signals 
derived  from  the  monopulse  difference  channels,  COLL  and  SKEW  correct  for 
collimation  and  non-perpendicularity  between  azimuth  and  elevation  axes 
and  ABIAS  corrects  for  misalignment  of  the  azimuth  reference  direction. 
These  corrections  wei-e  made  in  real-time  by  the  CAST  program  and  the 
required  auxiliary  data  and  model  parameters  were  recorded  for  further 
refined  post  mission  analysis.  The  raw  radar  data  together  with  the 
best  estimate  corrected,  edited,  and  smoothed  observations  used  for 
tracking  were  1’ecorded  for  post  mission  analysis.  The  smoothing  para- 
meters could  be  changed  easily.  For  the  observations  of  interest  to 
this  report  smoothing  was  accomplished  by  least  square  fitting  of  data 
from  successive  non-overlapping  20- second  periods  to  third  degree  poly- 
nomials and  reporting  the  value  from  the  fitted  curve  at  the  midpoint 
of  the  inverval.  L-band  tracking  data  obtained  during  the  joint  study 
and  the  first  half  of  the  1973-1974  low  elevation  angle  tracking  study 
were  obtained  using  an  eight-second  smoothing  interval. 

3.2  Orbit  Estimation 

The  smoothed  data  were  processed  by  the  multipass  orbit  determina- 
tion program,  ORBFIT  (Evans,  1973b),  to  determine  best  estimate  refer- 
ence orbits  for  use  in  estimating  refraction  effects.  Data  for  a single 
satellite  were  acquired  from  each  successive  pass  of  the  satellite 
visible  to  Millstone  Hill  radar  during  a five  to  seven  day  tracking  per- 
iod. The  data  were  from  different  directions  relative  to  the  radar 
station  and  represent  orbital  motions  in  different  directions  relative 
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to  the  radar  site.  The  smoothed  data  included  the  corrected,  observed 
position  estimate  for  each  smoothing  interval  and  the  rms  deviations  of 
the  raw  data  about  the  polynomial  curve  used  in  the  smoothing  operation. 
The  rms  data  were  used  to  establish  the  relative  quality  of  the  smoothed 
observations;  the  reciprocal  of  the  rms  deviation  value  was  used  to 
weight  the  observation  in  the  least  squared  determination  of  the  best 
estimate  state  vector  for  the  satellitr  orbit.  The  ORBFIT  program  also 
providod  best  estimate  fixed  bias  estimates  (EBIAS,  ABIAS)  as  well  as 
the  orbital  elements. 

ORBFIT  integrates  the  Keplarian  equations  of  motion  for  a satellite 
in  the  gravitational  force  field  of  the  Earth.  Although  capable  of  being 
considered,  ORBFIT,  as  operated  for  the  analysis  of  low  angle  tracking 
data,  ignored  the  effects  of  atmospheric  drag,  solar  radiation  pressure, 
and  perturbation  of  the  orbit  by  the  Sun  and  Moon.  The  Earth  and  gravi- 
tational force  model  employed  for  orbit  determination  was  the  DOD  WGS-66 
model  (DOD,  1966).  ORBFIT  used  differential  correction  to  estimate  the 
state  vector. 

The  use  of  ORBFIT  to  provide  precise  reference  orbits  was  verified 
by  comparison  with  the  Naval  Weapons  Laboratory  (NWL)  precise  orbital 
ephemeris  for  one  of  the  Navy  Navagation  System  Satellites  (NNSS) . The 
NNSS  were  used  as  the  beacon  sources  for  the  UHF  measurements  obtained 
during  the  joint  study  (Evans,  1973a).  The  NWL  data  were  derived  from 
multi-station  Doppler  tracking  data  and  were  processed  in  an  orbit 
determination  program  that  included  the  effects  of  drag,  solar  radiation 
pressure,  Sun-Moon  perturbations,  and  higher  order  geopotential  terms 
than  used  in  ORBFIT.  L-band  radar  observations  for  nine  passes  of  the 
NNSS  SPACETRACK  Object  No.  2965  from  the  8-11  August  1972  tracking  session 
were  processed  through  ORBFIT  for  comparison  with  the  NWL  data.  NWL 
ephemeris  data  at  one  minute  intervals  for  the  sections  of  the  orbit 
visible  to  Millstone  Hill  produced  the  following  rms  residuals  when 
compared  with  the  ORBFIT  reference  orbit:  range  92m,  azimuth  13  mdeg, 
elevation  1 mdeg.  When  ORBFIT  was  allowed  to  adjust  the  orbit  to  the 
NWL  data,  the  rms  residuals  were  reduced  to  20m,  1.6  mdeg,  and  0.3 
mdeg  respectively.  These  latter  numbers  are  a measure  of  the  ORBFIT 
program  analysis  errors;  the  former  are  a measure  of  the  state  of  radar 
calibration  in  1972. 


3.3  Calibration 


An  extensive  radar  calibration  program  was  undertaken  as  a part  of 
the  Lincoln  Laboratory-Bell  Telephone  Laboratory  joint  study  program  and 
was  continued  through  the  period  of  the  low  elevation  angle  tropospheric 
refraction  effects  study.  The  calibration  program  resulted  in  a contin- 
uous upgrading  of  the  error  models  for  the  facilities  and  in  the  obser- 
vation procedures  used  during  the  tracking  sessions . 

Angle-error  channel  calibration  data  were  obtained  on  each  satellite 
pass  that  rose  higher  than  about  10°  above  the  horizon.  After  the  first 
few  smoothed  data  samples  are  available  to  the  CAST  program,  an  optional 
ephemeris  update  can  be  performed  to  provide  best  estimate  position  data 
consistent  with  the  initially  entered  state  vector  for  the  orbit  and  the 
accumulated  observations  for  the  pass.  At  relatively  high  elevation 
angles  during  a satellite  pass,  error  signal  calibration  is  affected  by 
offsetting  the  antenna  in  elevation  then  in  azimuth  by  a predetermined 
angle  from  the  computed  best  estimate  satellite  location.  Offsets  are 
provided  in  both  positive  and  negative  directions  and  the  calibration 
values  for  the  error  channels  (EL  and  TRAV)  are  determined  from  the 
change  in  error  signal  versus  the  change  in  pointing  angle. 

Direct  azimuth  calibration  measurements  were  made  periodically 
using  optical  sightings  of  fixed  reference  targets  and  a boresight  tele- 
scope. The  azimuth  encoder  system  could  not  be  directly  mounted  on  a 
rotating  shaft  due  to  the  azimuth  rotaiy  joint  structure.  An  elaborate, 
octagon  within  an  octagon  cam  follower  system  was  constructed  to  provide 
a shaft  on  the  inner  octagon  that  rotated  once  per  revolution  of  the 
main  antenna.  The  accuracy  of  the  reported  azimuth  position  depended  on 
the  precisipn  of  the  location  of  the  rails  used  to  provide  the  outer 
octagon  surface.  This  |ystem  required  separate  modeling  of  the  errors 
introduced  by  the  eight  rails  in  the  system  requiring  the  determination 
and  monitoring  of  eight  separate  ABIAS  and  ASLOPE  calibration  constants. 

i r .y  ~ 

Analysis  of  a number  of  optical  target  sightings  revealed  that  the  cam 
follower,  encoder  system  had  an  rms  evrror  of  4.2  mdeg  about  the  calibra- 
tion model . 

The  tilt  meter  system  also  was  evaluated  by  direct  calibration  means. 
Each  meter  was  tilted  by  a known  amount  and  the  electronic  system  was 
adjusted  to  provide  a nominal  angle-to-volt  calibration  constant.  An 
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indirect  calibration  procedure  was  also  devised  to  measure  the  response 
of  the  tilt  meters  without  mechanically  disturbing  the  tilt  meters.  The 
antenna  was  rotated  while  pointed  in  a vertical  direction.  With  the 
antenna  in  this  orientation,  the  tower  should  not  be  deformed  by  antenna 
motion  and  the  output  of  the  tilt  meters  should  be  sinusoidal  functions 
of  azimuth.  The  departure  of  the  tilt  meter  outputs  from  the  expected 
sinusoids  was  taken  as  a measure  of  the  error  of  the  tilt  meter  system. 

In  particular,  d.c.  offsets  could  be  readily  detected  with  this  procedure 
as  well  as  changes  in  the  response  of  one  tiltmeter  relative  to  the 
other.  As  operated,  it  is  expected  that  tower  tilt  can  be  reliably 
measured  to  within  3 mdeg. 

The  overall  calibration  of  antenna  pointing  biases  were  made  using 
radio  star  observations.  Radio  star  positions  of  sources  suitable  for 
L-band  system  calibration  are  known  to  within  about  3 mdeg.  The  L-band 
monopulse  tracking  receiver  system  was  modified  to  enable  wide  bandwidth 
(1  MHz)  sum  and  difference  channel  radiometer  observations  of  the  radio 
stars.  The  computer  commanded  the  antenna  to  track  the  radio  sources  and 
to  provide  offset  calibration  scans.  Monopulse  tracks  of  the  radio 
sources  Cassiopeia  A,  Cygnus  A,  and  Taurus  were  periodically  obtained 
during  nighttime  hours  to  eliminate  interference  from  the  sun.  Calibra- 
tion measurements  made  with  sources  in  different  positions  in  the  sky 
were  combined  to  solve  for  the  coefficients  of  the  error  models:  for 
elevation  and  azimuth  biases,  collimation,  skew,  and  feed  droop.  Typical 
rms  uncertainties  in  the  determination  of  the  radio  star  positions  after 
removal  of  the  modeled  errors  are  the  order  of  5 mdeg  in  azimuth  and  3 
mdeg  in  elevation.  Radio  star  observations  were  made  periodically  dur- 
ing the  measurement  program  to  monitor  the  error  model  coefficients. 

Finally,  the  ORBFIT  program  provided  bias  error  estimates  as  well 
as  statistical  information  on  the  residual  differences  between  the  cor- 
rected and  smoothed  radar  data  and  the  estimated  orbit.  Some  of  the 
errors  such  as  feed  droop  could  be  monitored  by  repeated  runs  of  the 
program  with  different  values  for  the  parameter  to  be  investigated.  The 
best  value  for  the  parameter  would  consistently  yield  the  smallest 
residual  rms  deviations  for  several  different  satellites. 

Intercomparisons  between  the  ORBFIT  orbital  estimates  and  the 
smoothed  radar  observations  made  toward  the  end  of  the  joint  study  when 


the  calibration  procedures  were  most  advanced  and  when  the  best  models 
for  the  error  sources  were  used  yielded  rms  elevation  angle  residuals  of 
the  order  of  6 mdeg  and  azimuth  residuals  of  about  the  same  magnitude. 

In  particular,  for  all  the  satellites  tracked  during  the  27-30  March  1973 
tracking  session,  the  elevation  residuals  obtained  using  an  8-second 
smoothing  interval  were  6.1±0.3  mdeg  and  the  azimuth  residuals  were 
6.4±0.5  mdeg;  for  the  previous  tracking  session,  intercomparison  between 
the  ORBFIT  orbit  and  the  NWL  orbit  for  one  of  the  NNSS  yielded  rms  dif- 
ferences of  1.6  mdeg  in  elevation  and  3 mdeg  in  azimuth. 

During  the  1973-1974  time  period,  additional  modifications  were  made 
to  improve  system  calibration.  The  originally  provided  refraction  cor- 
rection coefficients  were  reintroduced  into  the  correction  procedure  to 
provide  a broader  based  statistical  correction  using  the  observed  sur- 
face index  values.  Tracking  with  fixed  offsets  was  introduced  to  provide 
higher  signal  to  noise  levels  for  the  calibrated  error  signals.  The 
operation  of  the  normalizing  amplifier  was  improved  to  reduce  uncertain- 
ties in  the  cross  section  value  and  in  the  error  signal  value.  Coherently 
detected  in-phase  and  quadrature  amplitude  values  were  also  recorded  for 
the  sum  channel  signal  to  provide  post  analysis  verification  of  the  nor- 
malization amplifier  AGC  output  and  to  provide  an  independent,  well-cali- 
brated, cross  section  measurement.  The  data  used  by  ORBFIT  was  con- 
strained to  lie  between  5 and  40°  elevation  angle  to  reduce  uncertainties 
caused  by  rapid  antenna  motion  at  higher  elevation  angles  and  by  refrac- 
tion effects  at  lower  angles.  This  last  step  removed  the  possibility  of 
having  the  best  fit  orbit  affect  the  determination  of  refraction  errors 
by  adjusting  the  orbit  to  compensate  for  any  errors  that  remained  after 

surface  index  correction.  The  final  results  for  tracks  of  0.1,  0.2,  and 

2 

1.0  m cross  section  spheres  were  rms  elevation  residuals  relative  to 

2 

the  reference  orbits  of  6. 0+0. 2 mdeg  for  the  1 m“  sphere,  6.8+0.2  mdeg  for 

2 2 
0.2  m spheres  and  7.4±0.2  mdeg  for  0.1  m cross  section  spheres. 

3.4  Sample  Observations 

Orbiting  spheres  were  used  to  provide  th.e  majority  of  satellite 
tracks  during  the  last  four  tracking  sessions  of  the  1973-1974  measure- 
ment period.  Sphere  data  were  selected  for  analysis  because  they  were 
not  subject  to  the  angular  errors  caused  by  glinting  (interference 
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between  scattered  signals  from  a complex  target  within  the  same  range, 
beam  width  resolution  cell).  Figures  1 and  2 provide  samples  of  the 
smoothed  output  data  provided  for  single  satellite  passes  by  the  CAST 
program.  Figure  1 is  for  a more  complex  target  with  a large  cross  sec- 
tion but  with  obvious  angle-of-arrival  fluctuations  at  all  elevation 
angles.  By  way  of  contrast,  Figure  2 displays  the  CAST  output  data  for 
a 1 m*  cross  section  sphere  (the  Lincoln  Calibration  Sphere,  LCS-4, 
SPACETRACK  Object  No.  5398).  The  cross-section  and  angle-of-arrival 
fluctuations  are  not  evident  except  at  low  elevation  angles.  Since  the 
sphere  does  not  adversely  affect  the  cross-section  or  angle-of-arrival 
measurements,  it  provides  an  ideal  target  for  refraction  effects  analysis. 
Unfortunately,  the  cross-sections  of  the  spheres  available  as  targets 
are  relatively  small  and  considerable  smoothing  (20  seconds)  is  required 
to  reduce  the  effects  of  receiver  noise. 

Two  types  of  data  vere  generated  by  the  CAST  program,  smoothed 

data  used  by  ORBPIT  for  orbit  determination  and  raw  data  for  use  in  post 

processing.  Figure  3 presents  a sample  output  for  a one  minute  segment 

2 

of  a track  of  a 1 m cross-section  sphere  at  low  elevation  angles.  The 
figure  displays  the  cross-section  values  measured  in  decibels  relative 
to  a 1 square  meter  cross-section  value  (dBsm)  and  the  elevation  and 
traverse  [azimuth  cos (elevation)]  angle  differences  between  the  refer- 
ence orbit  and  the  corrected  radar  return  versus  elapsed  time  and  true 
(corrected)  elevation  angle.  Data  for  each  radar  return  (each  pulse) 
are  displayed.  The  rapid  angle  fluctuations  are  caused  by  receiver  noise. 
The  expected  rms  angle  fluctuation  is  given  by 

A0  » 2 (!) 

n (S/N) 

where  0q  is  halfpower  beamwidth  (one  way),  n is  the  number  of  independent 
samples  averaged  to  determine  the  pointing  angle,  0,  (S/N)  is  the  signal- 
to-noise  ratio  and  A0  is  the  rms  angle  uncertainty.  For  the  data  dis- 
played on  the  figure,  the  rms  angle-of-arrival  fluctuations  are  between 
20  and  30  mdog  rms.  From  equation  (1),  the  signal-to-noise  ratio  for  a 
beamdwidth  of  0.7  degrees  and  a 30  mdeg  rms  noise  induced  fluctuation  is 
27  dB.  The  sum  channel  S/N  at  the  same  time  was  the  order  of  35  dB. 

With  the  tracking  offset  used  to  increase  the  signal  level  in  the  error 
channels,  a 27  dB  S/N  for  the  error  channel  is  reasonable.  Using  a 
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CAST  sumary  output  data  for  track  of  cylindrical  object  SPACETRACK  object  Ho.  727  obtained 
23  Septenber  1974  at  2048  UT 
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ELEVATION  ANGLE  ( DEG) 


20  second  smoothing  interval,  the  rms  angle-of-arrival  fluctuations  are 
reduced  by  a factor  of  4,5  from  30  to  6.7  mdeg. 

The  rms  angle-of-arrival  fluctuations  for  a typical  pass  of  the 

2 

1 m cross  section  sphere  are  presented  in  Figure  4.  The  apparent  ele- 
vation angle  to  the  sphere  is  listed  on  the  figure.  The  rms  values  are 
deviations  of  the  raw  data  from  a third  degree  polynomial  fitted  to  8.5 
seconds  of  data  (128  samples)  reported  for  the  midpoint  of  each  interval. 
Successive  data  points  are  for  overlapping  intervals  spaced  by  4.3 
seconds  (64  samples) . The  rms  values  at  angles  between  5°  and  40°  ele- 
vation angle  used  to  derive  the  reference  orbit  have  rms  fluctuations  of 
less  than  20  mdeg  rms.  The  20-second  samples  therefore  have  less  than  a 
4.5  mdeg  rms  fluctuation  produced  by  receiver  noise.  Over  much  of  a 
typical  arc,  the  rms  values  are  less  than  10  mdeg  producing  20-second 
averages  with  less  than  2 mdeg  rms  uncertainty.  Larger  elevation  angle 
fluctuations  are  apparen-  at  elevation  angles  above  20°  for  this  parti- 
cular pass.  These  are  sometimes  evident  when  the  antenna  moves  rapidly 
in  elevation  and  are  generally  suppressed  when  too  large  by  the  inverse 

rms  fluctuation  weighting  and  editing  employed  in  ORBFIT. 

2 

The  rms  angle-of-arrival  fluctuations  for  a 0.1  m cross  section 

target  are  displayed  in  Figure  5.  in  this  case,  the  rms  fluctuations 

2 

are  larger  than  for  the  1 m sphere  for  similar  observing  conditions 
(range,  elevation  angle).  For  this  particular  sample,  the  rms  fluctua- 
tions at  a 6°  elevation  angle  are  45  mdeg,  three  times  the  observed  value 

2 2 
for  the  1 m cross  section  sphere  (for  0.1  and  1 m spheres  at  the  same 

range,  the  ratio  should  be  3.2),  contributing  an  rms  uncertainty  of  10 

mdeg  to  the  measurement  of  the  20-second  average  angle-of-arrival  at  a 

2 

6°  elevation  angle.  At  10B  elevation  angle,  the  0.1  m cross  section 
sphere  would  contribute  6.7  mdeg  rms  fluctuation  to  the  20-second  aver- 
age angle-of-arrival  value. 

The  departure  of  the  elevation  angle  from  the  expected  (reference 
orbit)  values  at  low  elevation  angles  is  evident  in  Figure  3.  For 
elapsed  times  less  than  45  seconds,  the  elevation  angle  residuals  dis- 
play a relatively  smooth  trend  away  from  the  expected  value  with  the 
large  noise  induced  fluctuations  superimposed  on  the  trend.  At  five 
seconds,  the  average  deviation  is  8 mdeg;  at  40  seconds,  the  deviation 
is  24  mdeg.  The  smoothed  deviation  represents  the  uncorrected  component 
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Figure  4 RMS  angle-of-arrival  fluctuations  for  a high  elevation  angle 
pass  of  the  1 ra^  cross  section  sphere  SPACETRACK  object  No. 
5398  obtained  27  August  1974  at  1129  UT 
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of  the  bias  refraction  error,  the  refraction  induced  slow  departure  of 
the  surface  refractivity  corrected  elevation  angle  values  from  the  ref- 
erence orbit  positions.  The  20-second  smoothed  values  provide  estimates 
of  the  bias  refraction  effects  contaminated  by  the  residual  noise  fluc- 
tuations that  remain  after  averaging  and  residual  noise  fluctuations 
caused  by  tropospheric  angle-of-arrival  scintillation. 

The  data  displayed  in  Figure  3 indicate  that  the  satellite  set 
(disappeared  over  the  local  horizon)  at  about  53  seconds  elapsed  time. 

At  this  time,  the  cross  section  value  was  reduced  to  the  noise  level  and 
large  random  angle-of-arrival  values  are  evident.  In  the  five-second 
interval  prior  to  set,  large  amplitude  elevation  and  cross  section  fluc- 
tuations are  evident.  This  behavior  is  similar  to  the  results  reported 
by  Crane  (1976b)  for  the  low  elevation  angle  measurements  at  Haystack. 

The  more  extensive  instrumentation  used  for  the  Haystack  observations 
allowed  the  identification  of  these  large  amplitude  fluctuations  with 
internal  atmospheric  multipath,  a refraction  induced  phenomena,  not 
reflections  from  the  Earth's  surface.  The  relatively  large  amplitude 
fluctuations  together  with  the  elevation  angle  fluctuations  and  lack  of 
azimuth  angle  fluctuations  suggest  that  these  fluctuations  are  of  similar 
origin. 

2 

Cross  section  observations  for  the  1 m cross  section  sphere  are 
depicted  in  Figures  6 and  7.  In  these  figures,  the  values  are  plotted 
vs  apparent  elevation  angle.  The  cross  section  values  vary  about  the 
expected  cross  section  value  of  0.0  dBsm  reduced  slightly  by  the  effect 
of  atmospheric  focusing  (the  focusing  loss  is  expected  to  be  .6  dB  at  3°, 
see  Crane,  1976a;  see  also  next  section) . The  expected  cross  section 
values  including  focusing  loss  are  displayed  by  the  smooth  curves  on  each 
figure.  At  the  apparent  elevation  angle  values  listed  on  the  figure,  the 
one-way  gain  of  the  antenna  in  the  direction  of  the  surface  is  more  than 
18  dB  below  the  on-axis  gain  (Crane,  1971)  producing  a surface  reflection 
caused  variation  in  cross  section  of  less  than  ±2.5  dB  at  an  apparent 
elevation  angle  of  1,3°  and  of  less  than  ±0.5  dB  at  a 3°  apparent  eleva- 
tion angle.  The  observed  cross  section  fluctuations  are  significantly 
larger  being  more  than  10  dB  above  the  expected  value  at  an  apparent 
elevation  angle  of  1.06°  and  3 dB  above  the  expected  value  at  an  apparent 
elevation  angle  of  4°.  Refractivity  disturbed  conditions  are  clearly 
contributing  to  these  fluctuations. 
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Figure  7 Radar  cross  section  values  displayed  for  each  pulse,  rise  of 
1 cross  section  sphere  SPACETRACK  object  No.  5398  obtained 
27  August  1974  at  0953  UT 


The  Haystack  observations  suggest  that  a monopulse  tracker  would 
produce  elevation  angle  error  signals  that  were  a significant  fraction 
of  the  halfpower  beamwidth  of  the  antenna  in  response  to  internal  atmo- 
spheric multipath  events  producing  cross  section  increases  in  excess  of 
5-6  dB.  From  the  data  in  Figures  6 and  7 the  time  scales  of  these  fluc- 
tuations should  be  of  the  order  of  5 to  10  seconds;  the  effects  of  the 
fluctuations  should  be  reduced  to  variations  of  the  order  of  100  mdeg 
or  less  for  the  Millstone  radar  system  and  for  a 20-second  smoothing 
interval.  The  differences  between  two  smoothed  samples  obtained  from 
the  same  pass  of  the  1 m sphere  and  within  a 1°  true  elevation  angle 
range  were  as  large  as  60  mdeg  for  elevation  angles  between  1 and  2°, 

25  mdeg  for  elevation  angles  between  4 and  5°,  and  10  mdeg  for  elevation 
angles  between  9 and  10°.  These  changes  are  significantly  larger  than 
expected  due  to  signal- to-noise  limitations  and  indicate  that  the  fluc- 
tuations are  caused  by  either  internal  atmospheric  multipath  or  other 
small  scale,  large  gradient  refractive  phenomena. 
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4.  RAY-TRACING  ANALYSIS 


4.1  Ray-Tracing  Calculations 

The  coefficients  of  the  linear  relationship  between  surface  refrac- 
tivity, Ng,  and  elevation  angle  error  originally  provided  for  use  by  the 
joint  study  and  used  for  elevation  angle  correction  during  the  1973-1974 
measurement  period  were  derived  from  ray-tracing  calculations  using 
refractive  index  profiles  obtained  from  radiosonde  measurements.  The 
ray-tracing  program  used  for  the  analysis  was  developed  by  the  author 
while  at  Lincoln  Laboratory  during  the  1966-1968  time  period.  Since  some 
of  the  details  of  the  tracing  procedure  differ  from  ray-tracing  programs 
developed  by  others  and  the  program  has  not  been  documented  elsewhere, 
it  is  briefly  described  in  this  section. 

The  equations  used  for  ray- tracing  may  be  obtained  from  Maxwell's 
equations  (Freehafer,  1951),  or  from  variational  principles  (Bom  and 
Wolf,  1964).  From  either  starting  point,  the  basic  equation  of  geometri- 
cal optics,  the  eikonal  equation  may  be  obtained: 

VS  • VS  = n2  (2) 


where  S(^)  = path  function,  S(^)  = constant  being  the  equation  of  a wave- 
front,  n(£)  = index  of  refraction  (assumed  to  be  a real  (lossless) 
scalar,  n * 1 + N x 10~6,  N - radio  refractivity  or  refractivity) . The 
unit  normal  to  a wavefront  (unit  vector  along  the  ray)  is: 


s = 


VS 

n 


ds 


(3) 


where  £ = position  vector 

s = distance  along  the  ray. 


Then 


/ nds  = /VS  • d£  = S(£2)  - SQrp  = L 
ray  ray 


(4) 


where  L = path  length  between  points  ^ and  £2  °n  the  ray.  Equation  (4) 
provides  the  basic  equation  for  estimating  changes  in  transit  time  or 
electrical  path  length  due  to  the  atmosphere. 
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The  ekional  equation  may  be  obtained  from  equation  (4)  by  requiring 
that  the  path  length  be  stationary  along  the  ray  (Fermat's  Principle). 

The  equation  for  the  ray  path  may  be  obtained  by  finding  the  change  in 
7S  with  distance  along  the  path  using  equations  (2)  and  (3) : 

a?  <"  ■ It  <7S>  * § • 7(7S>  ■ ?7S  • 7(7S)  * 7" 

and  the  required  equation  is 

a?  (”  a?>  - 7"  <5> 

Equation  (5)  may  be  recast  as  a non-linear  coupled  set  of  first  order 
differential  equations  for  the  ray  path.  Using  the  following  geometrical 
relationships  to  describe  the  ray  trajectory: 


dr 

d?  = Sina 

(6a) 

d0  _ cosa  cos£ 
ds  r 

(6b) 

d#  m cosa  sin£ 
ds  r sin0 

(6c) 

where  a * local  elevation  angle  at  r,  0,  $ (angle  between  the  tangent 
to  the  ray  and  its  projection  in  a horizontal  plane  through 

£ * local  azimuth  angle  at  r,  0,  $ (orientation  of  the  projec- 
tion of  the  tangent  to  the  ray  in  the  horizontal  plane 
measured  clockwise  from  north), 

r,  0,  <p  * point  ;£  on  ray  trajectory  in  an  Earth  centered  spherical 
coordinate  system, 

equation  (5)  becomes 


da 

3? 


nr 


[cosa 


9 (nr) 
9r 


sina  (• 


cos£ 

r 


3 (nr)  + sin£  3 (nr). , 
90  r sin©  9$ 


(6d) 


d£ 

ds 


1 ,sin£  9 (nr)  _ cos£  3 (nr) . 
nr  cosa  *■  r 90  r sin©  9$  * 


cos<fr  cosa  sing 
r sin© 


(6e) 
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The  coupled  equations  for  the  ray  path  must  be  solved  numerically. 
For  special  cases,  the  equations  may  be  considerably  simplified.  For  a 
horizontally  stratified  atmosphere,  the  index  of  refraction  is  a function 
only  of  the  distance  from  the  center  of  the  earth.  Equations  (6a)  to 
(6e)  show  that  for  this  condition,  the  ray  trajectory  is  confined  to 
great  circle  planes.  A first  integral  of  equation  (6d)  may  be  performed 
yielding  Snell's  Law  for  a spherically  symmetric  medium, 

nr  cosat  = K = constant  (7) 


Equation  (7)  represents  a first  order  differential  equation  for  the  ray 
path,  cosa  being  related  to  the  differential  change  in  direction  with 
position  along  the  ray.  The  spherically  symmetric  equations  may  be  re- 
duced to  quadrature  expressions  for  the  path  length,  path  position,  and 
to  bending  along  the  ray, 


AL  (h2,hj)  = L(h2)  - L(hj) 


ro+h2 

Vhl 

rolh2 


A 0 (h2,hj)  = 0(h2)  - ©Chj)  = J 


ro+1h 


(nr)2  dr 
r\/(iir) 2 - K2  ' 

K dr 

W(nr)2'-  K2' 


(8) 

O) 


Aa(h2,h1)  = a(h2)  - aChj) 


(rnih2)n2 

f K d(nr) 

(r^hj)n,  (mr)* ^ 


(10) 


t (h2,hj)  = A0(h2,h1)  - Aa(h2,hJ 


(ID 


r 


e(h2>h1)  = a(hj)  - tan 


-1 


(r  + h,) 
o 1 


panAf0(h2,h1)]  'sinfA0(h2,h1)](rQ+h2) 


(13) 


where  L(h)  is  the  electrical  path  length  as  given  by  equation  (4) 
at  height  h, 

0(h)  is  the  central  angle  (earth  centered)  describing  the 
ray  position  as  shown  in  Fig.  8 at  height  h, 

a(h)  is  the  local  elevation  angle  at  height  h 

x(h2,h1)  is  the  bending  angle  or  the  angle  between  tangents  to 
the  ray  at  height  h2  and  height  hj. 


31 


e(h2,h.)  is  the  elevation  angle  error  or  the  angle  between  the 

ray  tangent  at  h.  and  the  actual  direction  to  h,  from  h,, 
and  1 21 

h is  the  height  of  the  ray  above  a reference  sphere  of 

radius  r . 
o 

The  quadrature  expressions  for  the  ray  may  be  explicitly  determined 
only  when  the  variation  of  n with  height  is  specified.  For  a linear 
change  of  n with  height,  equations  (8)  to  (10)  are  elliptic  integrals. 
Although  elliptic  integrals  have  been  tabulated,  the  changes  in  L,  0, 
and  o are  too  small  to  effectively  use  the  tables.  Alternative  approxi- 
mations also  may  be  used.  Explicit  evaluation  of  equations  (8)  to  (10) 
are  possible  when  the  product  (nr)  is  assumed  to  vary  linearly  with 
height.  For  m - nr,  Snell’s  Law  is  given  by  m cosa  » constant  and  the 
integrals  may  be  readily  evaluated  yielding: 


Vm2-K2'  + iv 


L , K , /?m  -K  + Vw  . v \ 

AMVhp-.-.cJ -) 


; w > o 


(13) 


K . -1  ( vr  + 2w  \ 
— T sm  ( — ;==========) 

y|w|  ' r/v2  - 4 uw' 


; w < 0 


where 

2 2 
w = b - K 

2 

u = a 
v = 2ab 

and  m=ar+b  or  n = a + — between  h.  and  h,; 

r it 
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j.hj)  = \fm 


AL  (h2,hj  = i/m2  - K2 


+ y-  A0(h2.hj)  + 


b log  (2ur  + v + 2a\/m2  - k2) 


(14) 


and 


Aa  (h2,h,)  - cos'1  (£) 


(15) 


Using  the  linear  change  in  m approximation,  measured  index  of  refraction 
profiles  (functions  of  height)  may  be  approximated  by  a number  of  layers 
each  with  linear  changes  in  m and  rays  are  traced  through  the  layers 
using  equations  (13)  to  (15).  This  process  is  significantly  faster  than 
numerically  integrating  the  original  ray  equations.  Index  of  refraction 
profiles  for  use  in  estimating  the  effects  of  atmospheric  refraction  are 
generally  obtained  from  radiosonde  balloon  flights  and  are  available  only 
as  tables  of  radio  refractivity  N = (n  - 1)  x 106  and  height.  The  vari- 
ation in  n between  the  reported  levels  is  not  provided.  The  linear  m 
approximation  allows  rapid  computation  of  the  ray  trajectory  given  the 
tabulated  values.  The  tracing  program  package  developed  at  Lincoln 
Laboratory  could  perform  either  the  rapid  linear  m profile  calculations 
or  the  slower  numerical  integrations.  Comparison  between  computations 
made  using  the  linear  m approximation  and  numerical  tracing  using  equa-^ 
tions  (6a)  to  (6e)  shows  the  errors  to  be  small  fractions  of  a percent, 
significantly  smaller  than  the  errors  due  to  an  imprecise  knowledge  of 
the  profile. 

The  ray  position  ©(h);  path  length,  L(h);  and  bending  t(h2,h1)  are 
parameters  that  describe  the  trajectory  of  a ray  in  space.  The  variation 
in  intensity  (magnitude  of  the  time  averaged  Poynting  vector)  may  also 
be  computed  using  geometric  optics.  Since  the  direction  of  the  ray,  s, 
is  the  direction  of  energy  propagation  or  of  the  Poynting  vector  in  a 
scalar  medium,  the  principle  of  conservation  of  energy  requires  that 
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where  I ■ intensity. 

Therefore,  the  intensity  is  proportional  to  the  area  of  a wavefront 
enclosed  in  a narrow  bundle  or  tube  of  rays.  For  the  case  of  a horizon- 
tally homogeneous  atmosphere,  the  focusing  factor  or  ratio  of  computed 
intensity  to  intensity  for  free  space  propagation  to  the  same  distance 
may  be  given  by  a quadrature  formula  and,  for  a horizontally  homogeneous 
atmosphere,  by  an  explicit  relationship  similar  to  those  for  AL,  A0,  and 
Aa: 

ICh,)  S(h,)cos2(a(h  )) 

T-fl-r  ■ F(h2)  ■ -j - 2 (16) 

fl  2J  t sin(a(hQ) ) sin  (A0(h2,hQ))  D(h2,hQ) 


where  D(h2,hQ)  » sin(a(h9)) 


A0(h2,hQ) 


rot)'2 

'l’ 

o o 


K3  dr 


r 2 „2T372 

r[m  - K ] 


i 
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If(h2)  is  the  intensity  for  free  space  propagation  and 

S is  distance  along  the  ray  from  hQ  to  h2;  hQ  is  the  initial 
ray  height. 

The  effect  of  extinction  by  atmospheric  gases  or  condensed  water 
(cloud  particles  or  hydrometeors)  along  the  ray  path  may  be  included  in 
the  evaluation  of  intensity  by  computing  the  line  integral  of  the  extinc- 
tion cross  section  per  unit  volume  along  the  ray: 


I = F(h2) 


Ife 


/ S(h2) 

Jo  Bd£ 


(17) 


where 


is  the  line  integral  of  g,  the  extinction  cross  section 
per  unit  volume  along  the  path. 


For  a horizontally  homogeneous  atmosphere,  the  line  integral  of  g is 
given  by:  S(h9)  r0+h2 

g(r)mdr 


S(hJ  r0+h 

f'Wj 


(13) 
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4.2  Simulation  of  Refraction  Effects 

The  refraction  effects  calculated  by  the  ray-tracing  program  are 
position  errors,  bending  values,  and  attenuation  due  to  focusing.  A 
set  of  sample  calculations  are  given  in  Table  4.  The  geometric  parameters 
h,  a,  G and  the  derived  parameters  L,  AL,  s,  r and  F are  listed.  The 
central  angle,  G,  is  tabulated  as  surface  distance  which  is  defined  as 
rQG.  The  starting  height  was  156  m,  the  height  of  the  phase  center  of 
the  Millstone  Hill  L-ba.nd  antenna  above  mean  sea  level.  The  25  km 
height  is  the  top  of  a typical  profile,  1000  km  the  height  of  the  orbit- 
ing spheres,  and  6371  km  is  one  earth  radius. 

The  ray  tracings  listed  in  Table  4 were  calculated  using  the  radio- 
sonde data  for  heights  below  30  km;  the  January,  45°N  latitude  U.S. 
Standard  Atmosphere  (ESSA,  1966)  for  heights  between  50  and  100  km;  and 
zero  radio  refractivity  above  (N  = 0;  n = 1).  The  profiles  were  smoothly 
connected  and  ionospheric  effects  were  ignored.  The  sample  calculations 
show  that  the  bending  values  are  within  99  percent  of  the  value  for  a 
1000  km  height  by  a height  of  25  km  at  an  apparent  elevation  angle  (a) 
of  1°  (5  mdeg  difference)  and  are  within  97  percent  of  the  value  for  a 
1000  km  height  by  a height  of  25  km  at  an  apparent  elevation  angle  of 
10°.  The  elevation  angle  errors  change  more  slowly;  the  differences 
between  elevation  angle  error  and  bending  (elevation  angle  error  at  an 
infinite  height)  are  4 percent  (19  mdeg)  at  a 1000  km  height  and  an 
apparent  elevation  angle  of  1°  and  2 percent  (4  mdeg)  at  the  same  height 
and  a 10°  elevation  angle.  The  range  errors  increase  rapidly  at  low 
heights  but  only  very  slowly  at  heights  above  100  km.  The  slow  change 
at  heights  above  the  region  for  which  N is  zero  is  caused  by  the  change 
in  curved  path  geometry,  not  the  change  in  refractive  index  along  the 
path.  Focusing  losses  also  continue  to  change  at  higher  altitudes  due 
to  the  change  in  curved  path  geometry.  For  radar  applications  (two  way 
propagation)  twice  the  focusing  loss  value  should  be  used. 

To  simulate  refraction  errors  for  the  Millstone  Hill  L-band  radar 
system,  a number  of  ray  tracings  were  performed  using  data  from  the 
nearest  radiosonde  station  with  a similar  climate.  Data  from  Albany, 

New  York  were  used  for  the  analysis.  Although  the  radiosonde  station  at 
Portland,  Maine  is  closer,  it  was  not  used  because  of  its  proximity  to 
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TABLE  4 


SAMPLE  RAY  TRACING 

Portland,  Maine  Radiosonde  Sounding,  17  February  1969,  00  UT,  Ng  = 304.2 


h 

a 

e 

L 

A L 

e 

T 

F 

Elevation 

Surface 

Range 

Elevation 

Elevation 

Focusing 

Height 

Angle 

Distance 

Range 

Error 

Error 

Bending 

Loss 

(km) 

(deg) 

(km) 

(km) 

(m) 

(mdeg) 

(mdejO 

(dB) 

0.156 

1.000 

10.0 

3.124 

277.34 

277.77 

54.45 

216.15 

369.60 

0.13 

25.0 

4.961 

489.63 

491.17 

63.33 

304.00 

442.08 

0.28 

100.0 

10.030 

1053.9 

1065.66 

65.01 

380.68 

446.95 

0.45 

1000.0 

30.176 

3294.0 

3643.77 

65.58 

427.57 

446.95 

0.55 

6371.2 

59.994 

6609.8 

10973.22 

65.66 

440.52 

446.95 

0.58 

0. 156 

2.000 

10.0 

3.572 

203.83 

204.26 

38.01 

154.05 

260.98 

0.07 

25.0 

S.  254 

398.15 

399.69 

45.70 

227.88 

326.45 

0.17 

100.0 

10.177 

946.12 

957.92 

46.90 

287.63 

331.10 

0.29 

1000.0 

30.221 

3175.0 

3524.69 

47.10 

313.69 

331.10 

0.35 

6371.2 

60.009 

6487.4 

10850.82 

47.16 

327.26 

331.10 

0.37 

0.1S6 

3.000 

10.0 

4.212 

156.83 

157.29 

28.55 

116.52 

196.65 

0.04 

25.0 

5.709 

329.49 

331.07 

35.16 

178.86 

253.98 

0.11 

100.0 

10.418 

853. S5 

865.41 

36.13 

227.53 

258.32 

0.20 

1000.0 

30.296 

3064.0 

3413.80 

36.22 

250.52 

258.32 

0.23 

6571.2 

60.035 

6370.9 

10734.41 

36.25 

255.84 

258.32 

0.24 

0.156 

S.000 

■ 

10.0 

5.809 

104.20 

104.77 

18.63 

76.28 

128.39 

0.02 

25.0 

6.968 

237.96 

239.73 

23.59 

122.14 

171.86 

0.06 

100.0 

11.152 

703.60 

715.77 

24.31 

157.29 

175.52 

0.10 

1000.0 

30.334 

2858.9 

3209.13 

24.34 

171.46 

175.52 

0.12 

6571.2 

60.115 

6148.3 

10512.23 

24.35 

174.28 

175.52 

0.12 

0.156 

10.000 

10.0 

10.425 

54.63 

55.  S7 

9.77 

39.64 

66.63 

0.01 

25.0 

11.107 

133.34 

135.90 

12.63 

65.61 

91.67 

0.02 

100.0 

14.093 

465.55 

479.62 

13.08 

85.76 

94.06 

0.03 

1000.0 

51.627 

2415.3 

2769.29 

13.09 

92.63 

94.06 

0.04 

6371.2 

60.491 

5625.0 

9993.39 

13.09 

93.67 

94.06 

0.04 

0.156 

20.000 

10.0 

20.209 

26.88 

28.65 

5.02 

19.45 

32.68 

0.00 

25.0 

20.558 

67.14 

71.72 

6.54 

32.56 

45.39 

0.01 

100.0 

22.261 

256.56 

277.16 

6.79 

42.89 

46.67 

0.01 

1000.0 

35.662 

1746.7 

2123.14 

6.79 

46.18 

46.47 

0.01 

6371.2 

61.966 

4671.7 

9073.09 

6.79 

46.56 

46.67 

0.01 

9 
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the  coast.  A sample  of  273  soundings  from  the  months  of  February  and 
August  for  the  years  1966,  1967  and  1968  were  used.  Every  00  and  12  UT 
sounding  from  each  month  was  traced,  but  only  273  soundings  reached  a 
height  of  25  km.  Scattergrams  of  bending  values  versus  reported  surface 
refractivity  are  depicted  in  Figures  9 and  10  for  apparent  elevation 
angles  of  1 and  10  degrees,  respectively. 

4.3  Surface  Correction  Model 

The  bending  values  to  a height  of  25  km  were  used  to  model  the 
elevation  angle  errors  to  be  expected  in  observing  satellites  at  altitudes 
above  1000  km.  As  evident  in  Table  4,  the  bending  values  at  a height 
of  25  km  are  a few  percent  lower  than  the  values  for  heights  above  100  km 
but  the  elevation  angle  errors  are  also  a few  percent  lower  than  the 
bending  values  for  heights  above  1000  km.  At  elevation  angles  above  5°, 
the  elevation  angle  error  at  a height  of  1000  km  is  within  1 mdeg  of  the 
bending  value  at  a height  of  25  km.  Since  the  expected  rms  fluctuations 
of  bending  (or  elevation  angle)  is  larger  than  2 mdeg  for  apparent  eleva- 
tion angles  less  than  5°  increasing  to  a value  of  12  mdeg  at  1°  when  sur- 
face correction  is  used,  a more  complex  height  dependent  elevation  angle 
error  correction  model  was  not  considered  and  the  easier  to  apply  bending 
correction  model  was  used  for  the  estimation  of  refraction  effects. 

A linear  regression  analysis  was  performed  on  the  bending-to-25-km- 
height  values  and  their  associated  surface  refractivity  values.  Ray 
tracing  calculations  were  performed  at  apparent  elevation  angles  of  0.1, 
0.2,  0.5,  1.0,  2.0,  3.0,  5.0,  10.0,  20.0  and  50.0  degrees  for  use  in  the 
regression  analysis.  The  result  is  a set  of  regression  coefficients  for 
each  apparent  elevation  angle.  Logarithmic  interpolation  was  used  to 
provide  coefficients  at  0.5°  apparent  elevation  steps  for  angles  below 
6°,  2°  steps  for  angles  between  6 and  20°  and  10°  steps  above  20°. 

Linear  interpolation  was  used  to  provide  the  coefficients  for  the  appa- 
rent elevation  angle  of  the  data  to  be  processed  by  CAST  or  ORBFIT. 

Table  5 presents  the  results  of  the  regression  analysis;  Table  6 provides 
the  interpolated  coefficients  used  in  the  computer  correction  of  eleva- 
tion angle  error. 
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Ng  - SURFACE  REFRACTIVITY  (N- UNITS) 

Figure  9 Simulated  elevation  bending  vs  surface  refractivity  for 

tracings  to  a height  of  25  km  and  an  apparent  elevation  angle 
of  1 degree,  Albany  data  set 
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Surface  refractivity:  Mean  = 326.6 
RMS  deviation  about  mean  = 21. 


TABLE  6 


5.  ANALYSIS  OF  REFRACTION  MEASUREMENTS 


5.1  Elevation  Angle  Error  Observations 

Radar  observations  of  spheres  obtained  during  the  last  two 
tracking  sessions,  9-14  September  and  23-28  September  1974,  were  analyzed 
for  refraction  effects.  These  data  represent  the  best  data  from  the 
1973-1974  measurement  period  having  both  the  most  accurate  calibration 
coefficients  and  the  most  refined  observing  techniques.  The  surface 
refractivity  values  for  each  tracking  session  are  displayed  in  Figures 
11  and  12.  The  refractivity  data  were  smoothed  for  1000  samples  before 
recording.  The  variations  of  surface  refractivity  over  the  eight  days 
of  observation  nearly  span  the  variation  encountered  in  the  Albany  data 
set.  The  diurnal  variation  in  Ng  is  clearly  evident  in  the  data  as  well 
as  longer  term  trends. 

The  smoothed,  corrected  radar  observations  were  compared  with 
the  ORBFIT  reference  orbits.  The  resulting  angle-of-arrival  differences, 
the  residual  elevation  angle  errors,  were  tabulated  in  1°  true  elevation 
angle  intervals.  Although  refraction  effects  may  vary  rapidly  over  ele- 
vation angle  intervals  as  small  as  1°  at  low  elevation  angles,  the  1° 
steps  were  used  for  analysis  since  the  major  elevation  angle  dependence 
had  been  removed  using  the  surface  refraction  correction  model.  The 
tabulated  values  for  all  sphere  passes  are  displayed  for  selected 
elevation  angles  in  Figures  13-15.  The  lines  displayed  on  each  figure 
are  the  results  of  a regression  analysis  of  residual  error  vs.  N . 

In  each  case  (see  Table  1)  the  correlation  coefficient  for  the  residual 
error  vs.  N was  not  significantly  different  from 

zero  (all  but  one  case  at  the  lo  level,  no  cases  at  a lower  significance 
level).  The  lack  of  correlation  implies  that  the  surface  correction 
procedure  using  the  apriori  simulation  results  removed  all  the 
variance  associated  with  Ng ; no  further  improvement  is  possible  using  a 
linear  relationship  between  elevation  angle  error  and  Ng . 

The  observed  residual  elevation  errors  are  not  correlated; 
therefore,  the  refraction  correction  procedure  can  be  worked  backwards 
to  provide  estimated  of  the  total,  uncorrected  refraction  error.  The 
resuits  of  this  operation  are  displayed  in  Figure  16  and  are  listed  in 
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■ SEQUENTIAL  DAY 
SEPTEMBER  1974 


370 


GMT  SEQUENTIAL  DAY 
SEPTEMBER  1974 


2-3°  ELEVATION 


• • « 
» 0 


290  300  310  320  330  340  350  360  370 

SURFACE  REFRACTIVITY  (N-UNITS) 

Figure  13  Residual  elevation  angle  errors,  2-3°  true  elevation  angle 
interval 
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RESIDUAL  ELEVATION  ANGLE  ERROR  (mdeg) 
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5-6°  ELEVATION  ANGLE 


290  300  310  320  330  340  350  360  370 

SURFACE  REFRACTIVITY  (N- UNITS) 


Figure  14  Residual  elevation  angle  errors,  5-6°  true  elevation  angle 
interval 
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Figure  15  Residual  elevation  angl 
interval 
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ELEVATION  ANGLE  ERROR  (MDEG) 


1000 


100 


APPARENT  ELEVATION  ANGLE  (DEG) 
I 2 3 5 10  20  50 


RADAR  OBSERVATIONS 
* NO  CORRECTION 
l f RESIDUAL  AFTER 
Vl  CORRECTION 


EXPECTED  N 
RAY  TRACIN' 


TRUE  ELEVATION  ANGLE  (DEG) 

Figure  16  Comparison  between  averaged  simulated  or  expected  and  averaged 
observed  refraction  effects 


Table  1.  The  elevation  angle  residuals  were  averaged  for  each  true 
elevation  angle  analysis  interval.  The  bending  correction  was  then 
removed  to  provide  estimates  of  the  mean  elevation  angle  error. 

These  are  plotted  as  radar  observations;  no  correction.  The  mean  bending 
values  calculated  at  a height  of  25  km  were  then  approximately  adjusted 
to  estimate  the  elevation  angle  error  at  1000  km  height.  This  adjustment 
was  accomplished  by  calculating  the  difference  between  bending  error 
and  elevation  angle  at  a height  of  1000  km  and  a surface  refractivity 
value  corresponding  to  the  mean  Ng  value  in  each  elevation  analysis  inter- 
val (ses  Table  3).  Adjustment  was  not  made  for  the  differences  between 
bonding  calculated  to  a height  of  25  km  and  bending  calculated  to  a height 
of  1000 ^km. 

The  results  of  the  comparison  betvaon  the  expected  (simulated 
adjusted  to  elevation  angle  error)  and  observed  elevation  angle  errors 
are  displayed  in  Figure  a6  and  Table  1 . The  differences  between  the  two 
refraction  effects  are  listed  together  with  error  bounds  to  represent 
the  effects  of  observed  fluctuations.  At  true  elevation  angles  below  2°, 
the  refraction  correction  was  underestimated,  above  2°  the  correction 
was  overestimated.  At  10°  elevation  angle,  the  difference  between 
bending  to  25  km  and  bending  to  1O00  km  height  should  be  the  order  of 
3 mdg  suggesting  that  although  near  perfect  agreement  is  shown  at  these 
elevation  angles,  the  residual  uncertainty  is  closer  to  3 mdeg.  This 
larger  value  is  reasonable  considering  the  residual  uncertainties  in 
system  calibration,  the  orbit  determination  program  and  the  neglect  of 
other  phenomena  such  as  ionospheric  bias  refraction.  The  latter  effect, 
although  relatively  small,  at  1295  MHz  (L-band)  still  can  contribute  as 
much  as  2 mdeg  increase  in  elevation  angle  error  in  the  elevation  angle 
range  of  interest  at  local  noon.  In  the  evening,  ionospheric  bias 
refraction  contributes  less  than  0.3  mdeg.  At  lower  elevation  angles, 
the  simulated  values  tend  to  be  too  large  in  the  3-5°  true  elevation 
angle  range  but  not  by  a significant  amount.  For  elevation  angles  below 
2°,  the  differences  between  the  measured  and  simulated  refraction  errors 
are  reduced  to  near  zero  when  the  additional  correction  for  bending 
above  a 25  km  height  are  included.  In  general,  the  observed  and  simula- 
ted mean  refraction  errors  agree  within  the  approximately  3 mdeg  measure- 
ment error  of  the  Millstone  L-band  radar  system. 
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5.2  Residual  Errors  After  Correction 


On  average,  the  residual  errors  after  correction  are  zero;  the  dif- 
ferences between  the  observed  and  expected  values  were  smaller  than  the 
overall  measurement  uncertainty  of  the  Millstone  L-band  radar.  The  same 
zero  average  residual  error  would  have  been  accomplished  if  the  average 
refraction  error  was  used  for  correction  rather  than  the  more  complex 
surface  correction  procedure.  The  test  of  the  surface  correction  pro- 
cedure is  the  measurement  of  the  reduction  in  variance  of  the  residual 
error  associated  with  the  use  of  the  correction  procedure.  Again,  since 

the  residual  errors  after  correction  were  not  correlated  with  N , the 

s 

reduction  in  variance  (or  in  rms  deviation  about  the  mean  or  correction 
estimate)  could  be  estimated.  The  residual  receiver  noise  induced  ele- 
vation angle  fluctuations,  however,  were  relatively  large  and  the  effects 
of  receiver  noise  were  modeled  prior  to  estimating  the  reduction  in 
variance. 

The  sample  observations  described  in  Section  3.4  showed  that  the 
receiver  noise  caused  angle-of-arrival  measurement  fluctuations  were 
sphere  size  (signal-to-noise)  dependent.  The  effect  of  receiver  noise 
plus  the  statistical  errors  in  orbit  fitting  were  reduced  by  processing 
the  data  by  sphere  cross  section  value  and  subtracting  the  variance 
estimated  for  receiver  noise  plus  orbit  fitting  uncertainty  from  the 
variance  of  the  observed  residuals.  Since  insufficient  statistical 
information  was  available  to  do  the  noise  subtraction  by  satellite, 
range,  cross  section,  and  elevation  angle,  a composite  method  was  used. 
The  variance  data  were  tabulated  by  elevation  angle  interval  and  satel- 
lite cross  section  value.  The  residual  fluctuations  about  the  reference 
orbit  at  elevation  angles  between  5°  and  40°  were  used  as  an  estimate  of 
the  composite  noise,  receiver  noise  and  orbit  fit  uncertainties,  and  the 
mean  square  residual  value  was  subtracted  from  the  tabulated  value  for 

each  elevation  angle  increment.  The  rms  noise  values  were  6.0  mdeg  for 

2 2 2 
the  1 m target,  6.8  mdeg  for  0.2  m targets  and  7.4  mdeg  for  the  0.1  m 

target.  The  resultant  rms  deviations  about  the  surface  corrected  values 
are  displayed  in  Figure  17.  The  data  are  displayed  by  sphere  cross 
section  and  a composite  value  for  all  the  sphere  measurements  is  provided 
which  was  calculated  by  summing  the  noise  corrected  variances.  The  dif- 
ferences between  the  variance  estimates  at  elevation  angles  above  8°  is 
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a better  measure  of  the  uncertainty  in  the  rms  value  estimate  than  the 
statistical  confidence  limits  calculated  from  the  number  of  degrees  of 
freedom  in  the  sample  population. 

The  composite  rms  residual  elevation  variations  about  the  surface 
model  correction  is  also  plotted  in  Figure  18  together  with  the  increase 
in  variance  contributed  by  removing  the  surface  correction  (radar  obser- 
vations - no  correction) . The  ratio  between  the  corrected  and  uncor- 
rected rms  values  is  a measure  of  the  effectiveness  of  the  surface  cor- 
rection procedure.  At  1°  and  above  7°  true  elevation  angle,  the  observed 
variation  of  elevation  angle  error  (no  correction)  is  identical  with  the 
value  expected  from  the  simulation  analysis.  The  observed  rms  fluctua- 
tions at  elevation  angles  between  1°  and  7°  true  elevation  angle  are  as 
much  as  40  percent  higher  than  expected  mirroring  the  influence  of  short 
term  small  scale  refractive  steering  and  internal  atmospheric  multipath 
discussed  in  Section  3.4. 

Although  all  the  variations  of  elevation  angle  error  associated  with 
surface  refractivity  variations  were  removed  by  the  surface  correction 
model,  the  model  did  not  significantly  reduce  the  variance.  The  rms 
fluctuations  about  the  correction  were  only  a factor  of  about  1.5  smaller 
than  the  uncorrected  rms  refraction  fluctuations  rather  than  the  factor 
of  6 expected  on  the  basis  of  the  simulation  analysis.  The  short  term 
fluctuations  which  were  not  modeled  and  cannot  be  sensed  by  surface 
based  techniques  significantly  affected  the  performance  of  the  refrac- 
tion correction  model . 

Summary  results  for  the  summer  (June  - August)  and  fall  (September  - 
November)  time  periods  from  the  Haystack  X-band  observations  (Crane, 

1976b)  are  also  plotted  on  Figure  18,  These  residual  fluctuations  are 
relative  to  a best  fit  third  degree  polynomial;  any  bias  effect  has  been 
removed.  They  provide  a measure  of  the  minimum  possible  rms  fluctuation. 
The  Haystack  observations  are  for  5 minute  and  1 hour  observation  periods. 
The  IDCSP  satellites  used  as  beacon  sources  for  the  Haystack  observations 
were  in  14  day  orbits;  the  apparent  elevation  angle  changed  by  approxi- 
mately 1 degree  per  hour.  The  1 hour  fluctuations  should  correspond  to 
the  fluctuations  to  be  expected  within  the  1 degree  analysis  interval 
used  for  the  L-band  radar  observations.  From  Figure  18,  it  is  evident 
that  the  rms  variations  about  the  surface  correction  model  are  of  the 
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same  order  of  magnitude  as  the  one  hour  fluctuations  of  the  September 
197S  Haystack  data.  These  fluctuations  are  caused  by  small  scale  tur- 
bulence or  by  bouyancy  waves  on  the  inversion  layer  capping  the  planetary 
boundary  layer. 

The  small  scale  nature  of  the  perturbations  is  evident  in  the  cross 
section  variations  (see  Section  3.4  and  Figures  3-5)  which  accompany 
the  angle-of-arrival  fluctuations.  Large  scale  horiiontal  inhomogeneities 
such  as  modeled  by  Gallop  and  Telford  (1975)  could  produce  angle-of-arrival 
errors  not  correlated  with  surface  refractivity  but  could  not  produce  the 
accompanying  cross  section  variations.  The  uncorrected  variance  repre- 
sented by  the  deviations  of  the  radar  observations  about  correction  are 
composed  of  contributions  by  both  large  and  small  scale  irregularities. 

The  effect  of  the  large  scale  irregularities  may  be  reduced  by  the 
addition  of  radiometric  observations  as  proposed  by  Gallop  and  Telford; 
the  small  scale  fluctuations  cannot  be  reduced  by  such  a technique.  It 
is  noted  that  the  large  scale  irregularities  will  affect  observations  of 
point  and  distributed  sources  (such  as  the  sun)  in  the  same  manner.  The 
small  scale  irregularities  will  affect  only  observations  of  point  targets 
such  as  the  radar  observations  of  satellites  or  observations  of  beacon 
sources.  This  difference  between  observations  of  point  and  distributed 
sources  is  responsible  for  the  difference  between  the  results  reported 
here  and  the  results  reported  by  Anway  (1963).  Anway  reported  that  the 
correlation  coefficient  between  elevation  angle  error  observations  and 
surface  refractivity  at  an  elevation  angle  of  8°  was  0.87  for  observations 
of  the  sun.  The  radar  observations  result  in  a correlation  coefficient 
of  only  0.61  (see  Table  3),  the  smaller  correlation  value  being  caused 
by  the  small  scale  irregularities. 
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6.  CONCLUSIONS 


The  results  of  the  low  elevation  angle  tropospheric  refraction 
effects  program  are  that  the  ray- tracing  model  simulations  of  refraction 
effects  are  in  agreement  with  observations  of  the  average  observed  values 
The  variance  of  the  observations  about  the  simulated  correction  value 
based  on  measured  surface  refractivity  was  significantly  larger  than 
expected.  The  observed  variations  about  the  correction  model  were  also 
larger  than  previously  reported  by  Anway  (1963).  The  departure  of  the 
measurements  reported  by  Anway  from  the  surface  correction  model  can  be 
attributed  to  the  effects  of  large  scale  horizontal  variations  in  refrac- 
tive index,  the  problem  described  by  Gallop  and  Telford  (1975).  The 
increase  in  variance  about  the  correction  model  observed  in  this  study 
was  comparable  to  the  variance  about  the  mean  value  of  elevation  error 
observed  within  an  hour  and  previously  reported  by  Crane  (1976b) , The 
source  of  this  increased  variability  relative  to  the  results  reported  by 
Anway  is  small  scale  refractive  index  gradients  which  affect  observations 
of  point  source  targets  but  which  do  not  affect  the  observations  of  dis- 
tributed sources  such  as  the  sun. 

The  effects  of  small  scale  refractive  gradient  induced  elevation 
angle  errors  cannot  be  reduced  by  regression  model  correction  schemes 
using  surface  refractivity,  surface  based  radiometry  or  a combination  of 
both  techniques.  At  this  time  these  fluctuations,  the  rms  residuals 
after  correction  displayed  in  Figure  18,  represent  the  uncorrectable 
limitation  to  low  elevation  angle  radar  observations  of  targets  at 
heights  in  excess  of  1000  km  (removal  of  the  effects  of  large  scale 
horizontal  inhomogenities  will  not  significantly  reduce  the  variances 
displayed  in  Figure  18;  the  effect  of  large  scale  irregularities  have 
already  been  eliminated  from  the  Haystack  data  by  the  processing  proce- 
dures used  to  prepare  those  data).  It  is  further  noted  that  these  resi- 
dual errors  cannot  be  removed  by  correction  using  ray  tracings  of  simul- 
taneously acquired  radiosonde  profile  data  since  the  radiosonde  data  do 
not  contain  sufficient  information  to  evaluate  the  effects  of  small 
scale  refractivity  gradients  at  the  time  of  a radar  measurement. 
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of 
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and  engineering  support  within  areas  o£  technical  competence 
u provided  to  ESV  Program  0Uic.es  (PCs)  and  other  ESP 
dements.  The  principal  technical,  mission  areas  are 
communications,  electromagnetic  guidance  and  control , sur- 
veillance o&  ground  and  aeAo  space  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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